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SUBORDINATION OF THE RESOLVENT FOR A SUM OF RANDOM 

MATRICES 

. . . , V. KARGIN 

^ ! 
O 

CN I Abstract 

This paper is about the relation of random matrix theory and the subordination 

C/^ ', phenomenon in complex analysis. Previously, subordination was discovered 

^ [ in free probability theory by Voiculescu (fT3\) and Biane ([ 1]). We find that 

approximate subordination holds for certain ensembles of random matrices. 

. This allows us to prove an improved local limit law for eigenvalues and a 

Dm ' 

p ^ I delocalization result for eigenvectors of these ensembles. 

^' 

1. Introduction 

The operator-valued analytic function Ga (z) = {A — zl)^^ is called the resolvent of 
^ ' operator A. Much of the modern approach to random matrices is based on the analysis of how 

OO I the resolvent of an A^-by- random matrix A at behaves when is modified, for example 

i when a column and a row is added. The results of this analysis are an important ingredient 

I in recent investigations of eigenvalues of Wigner matrices (||4l). In this paper we investigate 

^ I what happens with the resolvent of a matrix if an independent rotationally invariant random 

matrix is added to it. We find that the resolvent of the sum is (approximately) subordinated to 
the resolvent of the original matrix. 

Subordination is an important concept in the theory of functions of complex variables. If 
^ ' / (z) and g (z) are two functions analytic in the upper half-plane = {z : Imz > 0} , then 

■ f (z) is subordinated to g (z) if there exists an analytic function uj (z) : C+ — )• C+, such 

that f {z) = g (z + UJ (z)) (see ||9l). In this definition, / (z) and g (z) can be vector or even 
operator valued. 

An indication of the relation between subordination and random matrices gleamed when 
Voiculescu discovered ( ifTSl ) that the Stieltjes transform of the free convolution of probability 
measures is subordinated to the Stieltjes transforms of each of the original measures. Since 
free convolution is closely related to sums of large random matrices, this indicated that subor- 
dination should hold for random matrices. 

Voiculsescu proved his subordination result under some non-degeneracy conditions, which 
were later removed by Biane in HI. Moreover, Biane proved a significant strengthening of 
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this result by extending it to resolvents of sums of free operators. Let us recall the form in 
which the subordination result was obtained by Biane. 

Let ^ be a von Neumann operator algebra with the normal faithful trace r, and let be a 
von Neumann subalgebra of A. A conditional expectation t(- is a weakly continuous linear 
map A ^ B determined by the following properties: (i) t{1\B) = 1, and (ii) if 61, 62 G B, 
then T{biab2\B) = biT{a\B)b2. 

Biane's result (cf. Theorem 3.1 in [1]) says that if two self-adjoint operators A, B £ A are 
free, then the following identity holds for their resolvents: 



where t{\A) denote the conditional expectation on the subalgebra generated by operator A, 
and S{z) is a function analytic in the upper half-plane C"*". This function S (z) depends only 
on A and B, maps C"*" to C+ and satisfies the relation S{ir])/irj — ^ as r/ — ^ 00. 
In other words, t{Ga+b {z) \ A) is subordinated to Ga {z) ■ 

Since the concept of freeness is closely related to asymptotic behavior of large independent 
random matrices, it is natural to ask whether subordination results hold in the context of 
random matrices. 

More precisely, let A ^ and B^ be two A^-by- Hermitian matrices, and let Hn = An + 
UnB^U'^, where Un is a random A^-by-A^ matrix with the Haar distribution on the unitary 
group. The resolvent of is Ghm {z) = {Hn — z)^^ and the resolvents of matrices An 
and Bn defined similarly. 

It is known that for large A^, the behavior of operators An and UnBnU^ resembles 
the behavior of free operators. In particular, we can ask the question whether Ghj^, (z) is 
subordinated to Gaj^ (z) and Gbj^ (z) for sufficiently large A^. 

Since matrices An and UnBnU^ provide only an approximation to free operators, a more 
cautious version of this question is whether subordination holds in an approximate sense. In 
other words, we are interested in the existence of an analytic function Sb,n {z) : C+ — )• O^v 
such that 



where il^v is a region that approaches C+ as A^ — )• cx), and ||i2A,Af (^;)|| — as A^ — >• 00 
Here E denote the expectation with respect to the Haar measure on unitary matrices Un- (Of 
course, a similar question can be asked about subordination with respect to Gb^^- Since all 
considerations are the same, we will treat only the case of Ga^^ ■) 

It turns out that a suitable candidate for such a function is the following expression: 



t{Ga+b {z)\A) = Ga{z + S{z)) 



(1) 



EGhjv (-z) = Ga^ (z + Sb,n (z)) + Ra,n {z) 



(2) 




E/b^ (z) ^ 
EruHj^ {z) ' 



where 



(z) 
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is the Stieltjes transform of Hn and. 

fB^ (z) := N-'Tr{Br,—^). 

Hn - z 

{Sa,n (z) and /^^ (z) are defined similarly.) 

For this choice of 5*5^ AT (z) . it is known ( lITOl ) that the norm of the error term ||iiA,Af {z)\\ — )• 
as — oo provided that Imz > 770 > 0. Besides, it is clear that Sb,n (z) is analytic in 
C+ since Ktuhj^ (z) is the Stieltjes transform of a probability measure and therefore has no 
zeros in C+. Hence the main question is whether Sb,n (z) maps C+ to C+. Our first result is 
as follows. 

Theorem 1. Suppose A]\f = B^^ and Imz > 0. Then lm.SB,N {z) > and equality holds if 
and only if An and Bn are a multiple of identity matrix: An = Bn = cIn- 

Proof: Definitions imply the following identity: 

1 



Sa,n (z) + Sb,n (z) 



KniH [z) 

By assumption, An = Bn, and therefore, 

1 



Sa,n (z) = Sb,n (z) = -- { z 



2 V EmH{z) 

Let T] = Imz > 0. If A G M, then {X — ir])^^ £ Dr,, where Dr, is the closed disc with 
the diameter [0, i/r?] . Since Em/f {z) is a weighted average of (A — irf)"^ , we conclude that 
KruH {z) £ D^. Therefore, Im (Em/f [z))^'^ < —rj and Im z + (Kmn (z))"^ < 0. There- 
fore, Im {Sa,n (z)) = Im {Sb,n (z)) > 0. _ 

The case of equahty can occur only if Emn (z) G dDn which implies that eigenvalues 
of niH (z) are concentrated at a single point on the boundary of and this point does not 
depend on Un- Therefore Hn = 2cIn for a real c. This implies that An = Bn = cIn- D 

For the more general case of non-equal An and BN,'we. will prove the following result. 

Theorem 2. Assume that max{||^jv|| , ||i?jv||} ^ K and that z = E + ir] £ := 
{z: z> 0,\z\ < R}. Then for all N > cir]'^, 

mill {Im {Sa,n (z)) , Im {Sb,n (z))} > 

Nt]' 

where c\,C2 > depend only on K and R. 

Intuitively, this theorem says that for all sufficiently large N, the imaginary part of the 
subordination functions is almost non-negative. It is possible that in fact a stronger statement 
holds, namely that min{Im(5'^7v {z)) ,lm{SB,N (z))} > for all z G C+ and all suffi- 
ciently large N (which perhaps depend on Imz.) The author knows neither counterexamples 
to this stronger statement nor its proof. The proof of Thm. |2] will be given in the next section. 

Theorems [T] and [2] allow us to improve the estimate of the error term in Q. 
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Proposition 1. Let Ra,n {z) be defined as in (|2]). Assume that max{||j47v|| , ||i?Ar||} < K 
and z = E + 17] £ A/j = {z : r] > 0,\z\ < R} . Then, 

(i) if A = B, then \\Ra,n iz)\\ < c/ {Nrj^) , 

(ii) if A / B, then \\Ra,n {z)\\ < c/ [Nr]^) provided that Nif > ci, 
where the constants c and ci depend only on K and R. 

Note that by using this result we can estimate Ra,n {z) at points arbitrarily close to the 
real axis provided that is sufficiently large. This fact allows us to prove a local limit result 
for the convergence of eigenvalue counting measure. 

Let us introduce additional defintions needed to formulate this result. 

Let > . . . > A]y ^ denote eigenvalues of H^, and define the spectral measure 

of Hn as /ii^^ := Ylk=i ^ (h^v)- Define /i^^ and /x^Ar similarly to /ij^^, and assume 
that fiAif ~^ l^a and ^ibn ~^ /i/sasA^— >-oo, where and /i/j are probability measures and 
convergence is the weak convergence of measures. Let A/"^* (E) be the number of eigenvalues 
of Hn in the interval /* = [E - r]* , E + t]*] . What can be said about A/",,. (E) / {2r]*N) 
when A^ — oo? 

If 7]* is fixed, then it is known (HH) and (HIl) that the limit approaches fia ffl fJ-p (/*) , 
where ffl denotes free convolution, a non-Unear operation on probability measures introduced 
by Voiculescu in his studies of operator algebras. Local limit theorems address the question 
about what happens if rj* is not fixed but approaches when N ^ oo. 

We will formulate an answer in Theorem [3]below after giving some more definitions. 

Let fjLa and fif^ be two probability measures with bounded support, and let nia {z) and 
(z) be their Stieltjes transforms.The system of equations: 

m{z) = 7na{z + Spiz)) , (3) 
m{z) = mi3{z + Sa{z)) , 




= Sa{z) + 3(3 (z) . 



has a unique solution {m (z) , Sa (z) , Sf^ (z)) in the class of functions that are analytic in 
C+ = {z : Imz > 0} and have the following expansions at infinity: 

m{z) = -z-^ + 0{z-'^) , (4) 
Sa,b{z) = 0(1). 

The function m (z) , which we denote m^, is the Stiltjes transform of a probability measure 
which is in fact the free convolution of measures jia and /x^. It is known (|[T1) that Im^Q, (z) 
and ImSfs (z) are non-negative if Imz > 0. 

Definition 1. A pair of probability measures is called smooth at E if the limits of 

functions Sa (E) := lim^j^o Sa (E + irf) and Sp (E) := lim^^o {E + i-q) exist and if both 
ImSaiE) and ImS 13(E) are strictly positive. 
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Definition 2. A point E &M.is called generic with respect to measures {fia, fJ-p) if the follow- 
ing equation holds: 

k{E) := ^ ^ + ^ ^ -{E + SUE) + SUE)f + 0. (5) 

m'JyE ^ S^E)) mi^{E + Sa{E)) 

We will say that {^la, l^'ji) is well-behaved at E, if {^a,^p) is smooth at E and E is 
generic with respect to (^q, /i/j). It is known that if = Z^/? then {fia, fJ-fs) is smooth at E if 
and only if the free convolution ffl fifs has a positive density at E. Hence, smoothness is a 
strengthening of the requirement that fia H fJ-p is absolutely continuous at E. The concept of 
genericity is needed to guarantee that the solution of the system (|3]) at £' is stable with respect 
to small perturbation in the system. It comes from the requirement that the determinant of the 
Jacobian matrix of the system be different from at ii^. 

Next, let fii and fi2 be two probability measures on the real line, and let Fi {t) and F2 (t) 
be their cumulative distribution functions. The Levy distance between measures /xi and /Li2 is 
defined by the formula: 

d-L (/ii, 1^2) = supinf {s >0 : F2{x - s) - s < Fi{x) < F2 (x + s) + s} . 

X 

The Levy distance is in fact a metric on the space of probability measures and the convergence 
with respect to this metric is equivalent to the weak convergence of measures. 
Now we are ready to formulate the local limit theorem. 

Theorems. Letl ^ rj* ^ N^^^^ , and assume that (i) max {didjAN ^ lJ'a),dL{lJ'BM ^ I^ib)} ~^ 
0, (ii) the pair of probability measures (ponf^is) is well-behaved at E, (Hi) max{\\Ai\i\\ , ||i?7v||} < 
K for all N . Let denote the density of S l^p- Then, 

AT,* [E) 



2r]*N 

in probability. 



Pm (E) 



Here notation r/* » N^^/^ means that rf /N^^/^ — ;> 00 as iV — ;> 00. 

The theorem improves the local limit law in Q, where it was found that it holds for the 
window size if ~ (log A^)^^^^ . The expected optimal result is rj* ~ N^^^'^ with arbitrarily 
small positive e, which is similar to the case of classical Gaussian ensembles and the case of 
Wigner/sample covariance matrices. The proof of Thm. [3] will be given in Section [3] 

The fact that approximate subordination holds not only for Stieltjes transforms but also 
for resolvents of the matrices has an important consequence. Namely, this fact implies a 
delocalization result for eigenvectors of matrices H^. 

Let I wi^^ I denote an orthonormal basis of eigenvectors of and let aI^^ be the 

I J a=l 

corresponding eigenvectors. Let {cj}^^ be the standard basis and let ui^^ (i) denote the i-th 
component of vector ui^^ . 
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Theorem 4. Assume that (i) max. {dL{fj, a^j, fJ-a),dL{^BN' ^^|3)} ^ N^'^^^, (ii) the pair of 
probability measures (fia, f^p) i^ well-behaved in an interval I, and ( Hi) max { 1 1 \\ , 1 1 i?jv 1 1 } < 
K for all N . Then there exist constants C and c which depend only on (/i^, l^is) , /, and K, 
such that for all a with aI^^ G /, and alii = 1, . . . , N, it is true that 



P 



Note that ^^^Ii 



(^) ' 



> CiV"^/^ log ivj < exp 



1 and the theorem says that the size of individual coordi- 



2 

nates (i) is not greate than N^^/^ with high probability (up to a logarithmic factor). 
Hence, the result of this theorem can be interpreted as the claim that the eigenvectors of the 
matrix Hj\[ cannot be localized on the set of directions with cardinality significantly smaller 
than A^^/^. It is customary to express this by saying that the delocalization length of eigenvec- 
tors is at least N^/^. 

The assumption that {jj-a, jJ-js) is well-behaved is essential in this theorem. Otherwise, we 
can choose An diagonal and Bn = On, the zero matrix for each N. Then Hn = ^at is 
diagonal and clearly the eigenvectors are localized. In this case, it is easy to calculate from 
system ([3]l that Sfs {z) = and therefore the pair {fia, fJ-i^) is not smooth. 

The result is not optimal. The expected localization length is N, that is, eigenvalues are 
likely to be completely delocalized. 

The question about delocalization of eigenvectors occurs in physics for eigenvectors of 
random Schrodinger operators on Z*^. In particular, for d > 3 it is an open problem to show 
the complete delocalization for small disorder. Recently, there was some progress on delocal- 
ization of eigenvectors in simpler models, namely in the case of random Wigner matrices and 
in the case of random band matrices. In the former case, the complete delocalization has been 
established recently (see m for a review) and the method is similar to the method that is used 
in this paper. In the case of band matrices, it is expected that for matrices with the band width 
W greater than \/iV, the case of complete delocalization holds. What was actually shown in 
this case is that the delocalization length is greater than W^^'^^^ ([3|). The method is based 
on quantum diffusion and different from the method that is used in this paper. 

The rest of the paper is organized as follows. Section |2]devoted to the proof of Theorem |2] 
and Proposition [U Section |3] derives the local limit law of Theorem |3l and Section |4]is about 
delocaUzation of eigenvectors. Section |5]concludes. 



2. Approximate subordination 

The proof of Theorem |2] and Proposition [U will be given at the end of this section. First, 
we need to derive several preliminary results. 

Note that if Ra,n (z) were equal to zero, then we could easily compute Sb,n (z) as 

- (EGhj,)-^ + An - z. (6) 
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Hence, the following result is useful. 



Proposition 2. Let X be a random Hermitian matrix and let rj > 0. Then all eigenvalues of 
E {X — if]) belong to the half-plane M,^ := {z : Imz > r]} .In particular, 

all eigenvalues of — (EGh^)^^ + — z belong to C"*" = {z : Imz > 0} . 



the matrix 



Proof of Proposition |2l Every matrix {X (uj) — ir])~^ is normal and its eigenvalues are 
on the border of a disc D^j with diameter [0,i/rj] . (Here uj is an element of the probabihty 
space.) Hence, by Lemma[T]below, the eigenvalues of E (X — irj)"^ belong to the disc D^. It 



follows that eigenvalues of 



E (X - ir]) 



-1 



-1 



are in H^. □ 



Lemma 1. Let Aj, j = 1, . . . ,m, be a family of normal (finite-dimensional) operators. Sup- 
pose that the eigenvalues of all Aj are contained in a disc D C C, and let H = ^ Pj^j ^ 
convex combination of Aj. Then all eigenvalues of H are contained in D. 

Proof: By subtracting a multiple of the identity operator from all Aj, we can reduce the 
problem to the case when disc D has its center at 0. Assume that this is indeed the case. Let R 
be the radius of D. Since the operators are normal, their norms are equal to the maximum of 
the eigenvalue absolute values. Hence < R. Hence, \\H\\ < X^Pj ||^j|| < R- It follows 
that all eigenvalues of H have absolute value < i?. □ 

It is known (see lITOl or proof of Thm. 7 in Q) that the error term in ^ can be written as 
follows: 

Ra,N ■■ = ^^^{An - z)GAr,{z + SB,Niz))EAAj, (7) 

^ ^ -EA^^, (8) 



where 



We can re- write formula (O as: 

^Gh^ = Ga^ {z + Sb,n {z)) {I + ^Ar, (z)) , 



where ^a^, ■= i^N - z) EA^^. 



Hence, 

(EGh^)-^ = (/ + ^>A^r^ {An-z~ Sb,n (z)) 
If y^iv is defined as (/ + "^An)'^ ~ then we can further re- write this: 

Sb,n (z) = - (EGh^)'' + An-z + Ya^{An-z- Sb,n (z)) , (9) 

which resembles formula ^ except for a perturbation term. We are going to estimate the size 
of the perturbation and to establish the following results. 
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Proposition 3. Assume that ma,x{\\A]\f\\ , ||i?Ar||} < K and that z = E + i-q £ A 
{z: z>0, \z\ < R} . Then for all N > r]-^, 

\\Ya^ {An-z- Sb,n iz))\\ < c/ [Nt]') , 

where c > depends only on K and R. 



R 



(10) 



Proposition 4. Let A be an invertible diagonal matrix (not necessarily Hermitian), and U 
be a random Haar-distributed unitary matrix. Let B be an arbitrary matrix and assume that 
the matrix G (s) := {A + sUBU*)^^ is bounded with probability Ifor all s £ [0, 1]. Then 
EG (s) is diagonal for all s G [0, 1]. 

We also have the the following well-known result. (See Thm 6.3.2 on p. 365 in i5\). 

Lemma 2. Let X be a diagonal matrix and A be an arbitrary matrix. Then for every eigen- 
value \i of X -\- A, there exists an eigenvalue A, of X such that Xi — Xi < ||^|| . 

Assuming the validity of these results it is easy to prove Theorem [2| 

Proof of Theorem |2t Take the basis in which Aj^ is diagonal. From Propositions [2] and 

|4]we conclude that — (KGhn)^^ + An — z is diagonal and its eigenvalues belong to the 

halfplane C+ = {z : Imz > 0} . From Proposition [3l 

WYaj, {An-z- Sb,n {z))\\ < c/ [Nt]') . 
Hence, formula Q and Lemma |2] imply that 

lm{SB,N (z)) > 



Nrf 



The proof for Im {Sa,n {z)) is similar. □ 

For the proof of Prop. |3] we need several estimates which we formulate as Lemmas. 
In all of them we assume that max{||747v|| , ||-BAr||} < K and z = E + irj G Ar := 
{z: z> 0,\z\ < R}. 



Lemma 3. Let X be a random Hermitian matrix, \\X\\ < K. ThenlmK 



.-1 



> 



c?7, where c depends on K and R only. In particular, 
depends only on K and R. 



(KmHpf (z)) ^ < c' /t], where c' 



Proof: We can write 



ImE 



:Tr- 



1 



X 



a — IT] 



E 



1 rp V 



{X-aY + 7]^ 



Since {X - af+ifl < {{K + Rf + R^^ /, hence [{X - af + ifl^ > {{K + Rf + R^^ I 
and we conclude that 



E 



N {X- af + r/2 



> crj. 



□ 
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Lemma 4. For a numeric c > and every > 0, 

P{\'^Hi^ {z) -EniHj^ iz)\ > 6} < exp 

and 



2^4 



c6 7]' 

\\B 



N\ 



(11) 



2^4 



c5 Vj' 



rN^ 1 + 



\B 



N\ 



P{\fB^ {z)-KfB^ {z)\>5}<exp 

For the proof see Prop. 4 in Q. 
Lemma 5. 

P{\\Aa {z)\\ >e}< exp [-ce^r]^N^] , 
and W^Aa {z)\\ < c/ (A^?]^) where constants are positive and depend only on K and R. 



(12) 



Proof: For any Hermitian matrix X and > it is true that 
by using Lemma |4] we can write 



P{\\{mH^{z)-W.mH^{z))GHA\>5h} < exp 
P{\\GA^{fB^{z)-EfB^{z))GHj>S/r]^} < exp 



{X — irj) ^ <l/?y. Hence, 



2:„4 



\B 



N\ 



B 



N\ 



Set e = d/f] and e = in the first and the second inequalities, respectively, and use the 
triangle inequalit for norms and we wiU get that 



P{\\AAiz)\\>e} < exp 



min { rf,rf/ [ 1 + 



\\Bn\\ 
< exp [-ce^T]^N^] , 



V 



B 



N\ 



where we used the assumption that \z\ < R and where the constant c in the last inequality can 
depend on \\Bn\\ and R. For the second claim, note that ||EAa|| < E || A^H by the convexity 
of norm, and E || A^H can be estimated by using the first claim of Lemma |5] and the equality 



/■oo 

EX= / il-Txit))dt, 
Jo 



valid for every positive random variable X and its cumulative distribution function Tx (t) ■ 
In our case, we obtain 



poo 

E ||Aa|| < / exp [-ct'^Tj^N^] dt 
Jo 



,4- 



n 



Lemma 6. ||^Ajv II ^ c/ (^Ntj^) , where c > depends only on K and R. 

Proof: This is a direct consequence of Lemmas [3] and [5] which estimate (Em/^^ 
||EA^^|| as c/rj and c/ (^Nr]^^ , respectively. □ 



and 
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Lemma7. If \\X\\ <e <l/2,then {I + X) 



\-i 



< 2e. In particular, for all N ^ rj 



II^Ajv II — [Nt]^^ ,where c > depends only on K and R. 

Proof: (/ + X)-^-I = X ^Zo ("1)' and || -X EZo ("1)' < ^ EZo < 
2e. The second claim follows from Lemma|6l □ 

Lemma 8. \SB,Niz)\ <c/rf'. 

Proof: By definition, Sb,n {z) = —K/sf^ {z) /Em/Zj^ (z) . From Lemma|3l (Em^jv (z))^^ 
c/rj. In addition, 

1 



Wb^ {z) 



E— TV f Bn ^ 

\ Hn-z 



< \\Bn\\^ 



Hn - z 



1 

< C-. 



□ 



Proof of Prop.|3j Estimates in Lemmas |7] and [8] imply (ITOl ). □ 



< 



Proof of Prop. |H For proof we need several Lemmas. 

Let Uij and u^j denote the matrix elements of random matrices U and U* , respectively, 
where U is distributed uniformly on the unitary group U {N) . Consider the product 



n 



Let L (i) be the number of times index i appears as the first index of Ui^j^ , and let R {i) be the 
number of times it appears as the second index in u^j^ . 

Lemma 9. IfL (i) ^ R{i) , then EH = 0. 

Proof: Let (i) be a diagonal matrix that coincides with the identity matrix everywhere 
except at the i-th place of the main diagonal where it has exp (t\^—lj instead of 1. Assume 
that t is random and distributed uniformly on [0, 27r] . For every t, AtU is distributed uniformly 
on U {N) , hence if Ht denotes the product analogous to 11 but computed for the matrix AfU, 
then the conditional expectation E {Ilt\t) = E (11) . Hence E (lit) = E (11) , where the first 
expectation is taken both over U and t. However, if L (i) / R (i) , then we have lit = 
gfetv^jj ^jjj^ ^ _^ Q_ jjgj^^g j^27r jj^^^ = and E (Ut) =0.0 

Lemma 10. Let A is a diagonal matrix (not necessarily Hermitian), and B is arbitrary. Let 
U be a random Haar-distributed unitary matrix. Let n be a non-negative integer Then 

E[{UBU*A)''UBU*] 

is diagonal. 

Proof: A typical matrix element of {UBU*A)"- UBU* is 
We need to show that the expectation of this expression is zero if ii / In+i- 
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By assumption, A is diagonal, hence every non-zero product has to satisfy the equaUty 
Is = is+1- Let ii = i. Then L{i) > R (i) since if i appears as the second index of tt^^/^, 
s < n, then it must appear as the first index of Ui^^-^j^^-^ , and it also appears as the first 
index of ui-^j-^ . Moreover, the equality holds if and only if i appears as the second index of 
^fc„+i«„+i- Lemma|9j the expectation of the product is zero unless L{i) = R (i) , that is, 
unless ii = In+i- D 

Now we can complete the proof of Prop. |4l Let 

G{s) = {A + sUBU*y\ 
If s is small, then we can expand G (s) as a convergent power series in s, 

G (s) = A-^ + A-^ If^ [UBWA-^Y 

\n=0 / 

By using Lemma [TOl we find that if s is sufficiently small and the series are convergent, then 
EG (s) is diagonal, that is, that all off-diagonal elements of EG (s) are zero. By assumption 
the matrix G (s) := {A + sUBU*y^ is bounded with probability 1 for all s G [0, 1] . Hence, 
G (s) can be analytically continued along this interval, and EG (s) is an analytic function for 
all s S [0, 1] . Off-diagonal matrix elements of EG (s) are zero for all sufficiently small s, 
therefore they are zero for all s e [0, 1] by properties of complex-analytic functions. Hence, 
EG (s) is diagonal for all s G [0, 1] . □ 

Proof of Proposition [H If A^ = Bj^, then Theorem [T] implies that 

Imz' = Im (z + Sb,n (z)) > r], 

where rj := Imz. Similarly, if A ^ B and N > clrf for sufficiently large c, then Theorem|2] 
implies that 

Im/ = Im (z + Sb,n (^)) > ??/2. 

Hence, in both cases ||Gajv (-^OII ^ c/r/. Since Ra,n {z) = Ga^ i^') (z) , we can 
use Lemma |6] in order to obtain 

\\RA,Niz)\\ < \\Ga^ {z')\\ W^A^ iz)\\ < 

□ 



3. Local Law for Eigenvalues 

We will prove Theorem |3] at the end of this section after several preliminary results. 
First, we can exclude mg (z) from system ([3]): 

ma{z + Sf^{z)) + , g, . ^ , g, . X = 0, (13) 

Z + Sa{z) + Si3{z) 
mp{z + Sa{z))+ , g, / X , g, / X = 0. 
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A similar system can be written for Sa,n (z) and Sb,n (z)' 

1 

z + Sa,n (z) + Sb,n {z 



ruAt^ {z + Sb,n (z)) + , , c. , c_ = i"A,N {z) , (14) 



me^ {z + Sa,n (z)) + -—^ ^ -— = vb^n (z) , 

Z + SA,N [Z) + bB,N (Z) 

where ta^n (z) = N-'^Ti {Ra,n (z)) , rB,N (z) = iV^^Tr {Rb,n (z)) , and Ra,n, Rb,n are 
defined as in (|7]). 

There is also an "intermediate" system which uses and m/3 as in ([T3] ) and has non-zero 
right hand side as in ([T4l) : 

nia (z + (^)) H ^ = rA,N (z) , (15) 

(z + Sa (z)) + ^— ^— = rB,N (z) ■ 



We are going to show that the solutions of (fTSl ) (^Sq (z) , S"^ (z) j are close to solutions of 
both ([13]) (5o (z) , Sp {z)) and {Sa,n {z) , S'b^at (z)). In this way, we will show that the 
solutions of ([T3] ) and ([T4l ) are close to each other. 

Proposition 5. Assume that (i) a pair of probability measures {^a, f^/s) well-behaved at E, 
(ii) max{||AAr|| , ||i?Ar||} < K, (Hi) z = E + ii] e {z : i] > 0,\z\ < R} , and (iv) Nr]^ > ci, 
where c\ depends on K and R. Then for some positive sq, tjo and c, ifdLilJ'A, l^a) < s < sq 
and dilfiB, fJ-p) < s < sq, then 



and 



for all rj < rjQ. 



\So, {E + ir]) - Sa,n {E + < + cs, 



\Si3{E + ii]) - Sb,n {E + ir/)| < + cs. 



Proof: The claim is the direct consequence of Lemmas [TT] and [12] below, in which the 
differences between solutions of (TSh and ( [T5] ). and between solutions of $T5\i and ([14] ) are 
estimated. □ 

In all lemmas below we assume that a pair of probability measures {^a,fJ-i3) is well- 
behaved at E, that max{||^Ar|| , ||i?Af||} < K, and that z = E + irj e {z : rj > 0,\z\ < R} , 

Lemma 11. Assume that Nrj^ > ci, where ci is a constant that may depend on K and R and 
the pair of measures (fia, fJ-js) ■ Then for some positive rjQ and c, it is true that 

c 



Sa {E + irj) -Sa{E + il]) 



< 



and 

Sfi{E + iri)-Sp{E + ir,) < 



for all T] < rjQ. 
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Proof: Let F (t) be a function defined by the formula: 

{z + Si + S2)~^ + ma{z + S2) - TA^N {z) 

{z + si + S2y^ + mp{z + Si) - rB,N (z) 

1/2 




F 



Let us use the norm ||(xi, X2)|| = + \X2\ 

The system ([15] ) can be written as F (x) = 0. By Proposition [T], \\F {Sa, Sp) 
c/ [Nrf') provided that Nrf > ci. 

The derivative of F with respect to s is 

/ _ _L _L 



{rA,N,rB,N)\\ < 



[Z + Sl + S2) 

-2 



{z + si + S2) ^ + rn'a {z + S2) 



{z + si + 82)' 



\ - {z + si + S2) + m'^ {z + si) 
The determinant of this matrix is 

[m'„ {z + S2) + m^Jj (z + si)] {z + si + ^2)"^ - m'^ {z + S2) m'^ {z + si) . 

This expression is non-zero at {Sa,Sis) by assumption ([5]). This, together with another ap- 
plication of the assumption of smoothness, shows that the entries of the matrix [-F']^^ are 
bounded at {Sa (z) , Sfs (z)) for all z in a neighborhood of E. This shows that operator norm 
of [F']"^ is bounded. 

Similarly, the assumption of smoothness implies that the operator norm of F" is bounded 
for all (si, S2) in a neighborhood of {Sa (z) , S"^ (z)) and all values of parameter z in a neigh- 
borhood of E. 

It follows by the Newton-Kantorovich theorem (||6l) that the solution of the equation 
F (s) = 0, which is (^Sa (z) , Sjj (z)^ exists provided that Nrf > ci and satisfies the in- 
equalities: 



Saiz)-Sa{z) <c/{Nrf) 



and 



□ 



Sr iz) - Sr (z) 



< c/ [Nr]^) 



Corollary 1. If N :^ r]^^ and r? < 1, then Sa (E) := lim^^o Sa {E + irj) and Sp (E) r- 
lim^j^o {E + ir]) exist and both lra.Sa{E) and ImS'^(£') are strictly positive. Moreover, 
1 1 



+ 



{E + Sa{E) + Sp{E)Y ^Q. 



m'^{E + Sp{E)) m'p{E + Sa{E)) 

Lemma 12. Assume that Nrf > c\ , where c\ is a constant that may depend on K and R and 
the pair of measures, (fia, f^/s) ■ Then for some positive sq, r/g and c, if dL[jiA-, l-ta) < s < sq 
and dL{^Bi l^i-p) < s < sq, then 



and 



Sa (E + irj) - Sa,n (E + ir]) 
Sfi {E + irj) - Sb,n {E + iri) 



< as, 



< cs, 
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for all r] < TjQ. 



Proof: We proceed similarly to the proof of Lemma[TT]and write system (O as F/v (x) 
0, where 

I + ^1 + •^2)"^ + "T'Ajv + S2) - ^A.Af (2) \ 
52/ V {z + Si+ S2)~^ +mB,N {z + Si) - TB^N {z) j 



F : 



Corollary [T] implies that [-F']^ is not at \^Sa, S'/jj . It is also possible to check that F" 

is bounded in a neighborhood of (^Sa, Sjs^ . Then we use the Newton-Kantorovich theorem 

to find a solution of ([141 ) by starting with the initial guess (^Sa, Sjs^ . In order to estimate the 
initial error, we can use the following lemmas. 

Lemma 13. let mi (z) and m2 (z) denote the Stieltjes transforms of measures ni and fi2, 
respectively. Let (/ii, ^2) = s cind z = E + ir], where ij > O.Then, 

(a) \mi (z) — 7712 (-2)1 < csrj^^ max |l, r/~^} where c > is a numeric constant, and 

(b) I ^ (mi (z) — 7772 iz))\ < Crsr]~^~'^ max {l, r/^^} where c > are numeric constants. 

This lemma was proved as Lemma 2.2 in |[8l . 

Lemma 14. Suppose that {fJ-a, l^A^) < ^ '^^^ {jip^fiBf^) < s. Then, 

1 



Z + Sa (z) + SfS (z) 



and 



+ ruA^ [z + Sjj{z)) - rA,N (z) 



+ mBff {z + Sa{z)) - rB,N {z) 



< cs, 



< cs, 



Z + Sa (z) + Si3 {z) 

for all T] > cN^^^^, where c is a positive constant that depends only on the pair of measures 
Ha and fijs. 

Proof: Sa{z) and S/^ (z) satisfy the equations of system ([TSl l. which implies that it is 
enough to show that 



rriAis, (z + Sfs {z)) -ma[z + Sp {z) 



and 



m-Bjv [z + Sa{z)] - mpiz + Sa {z) 



< cs, 



< cs, 



This follows from Corollary [T]and Lemma[T3l Indeed, let 770 (E) > denote m.m{lin.Sa (E) , ImS^ (E)}. 
Then, by Lemma [T3l 



"T-Ajv yz + Sp {z)j -ma\^z + Sp (z) 
and similar estimate holds for the difference 



< csmm 



7775^ [z + Sa{z)\ -mii[z + Sa (z) 



. □ 
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By using this Lemma and the Newton-Kantorovich theorem we find that the distance be- 
tween [Sa,Sj3] and the solution of (fT4l ) is smaller than cs. This completes the proof of 



|Em/f^ {E + IT]) -mm{E + ir])\ < + cs, 



Proposition 15] □. 

Proposition 6. Let the assumptions of Prop, \5\hold. Then for some positive sq, tiq and c, if 

dii^J-A, iJ-a) < s < So and dLi^J-B, ^J-fi) < s < sq, then 

c 

Nrf 
for all rj < r]Q. 

Proof: Since Kmnj^ = {z + Sa,n + Sb,n)^^ and = (z + Sa + Sjs)^^ , therefore 
Emj^^ - ma = — — — ^— — , „^ , o (^<^ + 3/3- Sa,n - Sb,n) ■ 

Z + JA,N + ^B,N Z + ba + 0/3 

The denominator is bounded away from zero if Nt]^ > ci by Prop. [5] (luiSa (E) and 
ImSfs [E) are bounded away from by the assumption of Prop. |5] and Sa,n and Sb,n are 
close to ImSa {E) and ImS*/? {E) , respectively, by its conclusion). The numerator can also 
be estimated by Proposition [5] □ 

Proposition 7. Suppose that the assumptions of Proposition\5\hold. Let rf = Mr] and Ir^* = 
[E — T]* + ir], E + r]* + ir/] . Then for some positve c and ci, that can depend on K, R, M, 



sup \mH,N (z) - mm iz)\ > + cs > < exp ( -ci (logiV) ) , (16) 



ifV < Vo- 



Proof: This is a combination of statements in Proposition [6] and in the following strength- 
ening of Lemma m which was proved as Corollary 6 in |7]. 

Lemma 15. For some positive c and c\ which can depend on M and for all (5 > 0, 

P < sup Im^jv {z) - EmHf^ {z)\ > 6> < exp ( - ^ iV^ 

provided that N > ci ^-y/— log {r]6)j / {rfS) ■ 

Let us take 5 = clog N/ {^Nif ^ . Then the assumption of Lemma [T5] is satisfied provided 
that 7] = O (A^~^+^) . In particular, if r/ = cN^^/^ for a sufficiently large c, then Lemma [T5] 
implies that 

pi sup \mH{z) - ¥.mH{z)\ > \ < exp (^-c(logiV)^ 



In addition, Proposition |6] says that 

I Em// 

Together, these statements imply the claim of Proposition |7] □ 



Em// (z) - msi {z)\ < + cs. 
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Proof of Theorem |3j The proof is similar to the proof of Corollary 4.2 (local semicircle 
law) in 0. We provide details for the convenience of the reader. Define r] = A^~^/^, and 
r]* = Mr]. Let 

R(X) : = - / ^ dx 

^JE^r [x-\f + rf 

^ ( (E-\ \ (E-\ ,^ 
= — arctan V M \ — arctan I M 

vr V V ^ / V ^ 

Then i? = 1/* + Ti + r2 + Ts, where Ti, Ts satisfy the following properties: 

|Ti| < c/^/m, supp(Ti) c Ji = - 27?*,^ + 2r?*] , 

\T2\ < 1, supp(r2) C Ji U J2, 
where Ji and J2 are intervals of length y/Mrj with midpoints at E — r]* and E + r]*, respec- 
tively, 

iTsI < ^, supp(r3) G II 

Note that 

^ - ^/lHA).„„(.A) 

= 2^/ / m + r2 + r3)M^^^(dA). 

The last integral can be estimated as O (M~^/^) . For the main term we have 

1 /• 1 1 f I 

/ — Immffl (x + ir/) dx H / —lm(mH,^{x + irj)—ms(x + irj))dx. 

2r]* J J* vr 2ry* Jj* vr 

The first part converges to (E) because the assumption that (/i^, /i^) is well-behaved at E 
implies that /i^ ffl /i^ has an analytic density in a neighborhood of E. For the second term we 
can use the estimate in Proposition |7] that shows that this term converges to in probabihty as 

N ^ oo.D 

4. Delocalization 

Proof of Theorem |4l For simplicity of notation, let us omit the subscript in Ghj^ , 
i.e., G := Ghj^- We can write 

Let us set E = Xa, then 

\G^^ {E + i7j)\> ImG,, (E + i7]) > ^^^iilL, 

V 

and therefore 

\va <rj\Gii {E + ir])\. 
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By using Lemma [TSl below with rj = N^^^^ log we obtain the claim of the theorem. □ 

Lemma 16. Let 1 S> » N^^/^, and assume that (i) max {didi a^^, lJ'a),dLil^Bff, IJ'/b)} ^ 
A^^^/^, (ii) the pair of probability measures (fia^fJ-i^) is well-behaved at E, and (Hi) max {\\ A , \\Bj\f\ 
K for all N . Then, for all k = 1, . . . ,N, 

\EGkk {E + iv)\ <C, 

where C is a constant that depends only on K, E, and the pair (^a, fip) . 

Proof: Recall that (as in © 

EGHr, (z) = Gam + Sb,n (z)) + Ra,n (z) , (17) 
By Proposition |5l assumptions of the lemma imply that 

IS^iE + irf) - Sb,n {E + ir])\ < + cs, 

where s = max{dL(^Ajv, Aio), c?L(/^Biv; M/?)} • Hence z — Sb,n [z] has a positive imaginary 
part for sufficiently large A^. In addition, by Proposition [H 

\\Ra,n (z)\\ <c/{Nrf). 

Hence, 

\\EGHAE + ir^)\\<C, 
which implies the claim of the lemma. □ 

Lemma 17. Let z = E + irj where t] > 0. Then, for a numeric c > and every 5 > 0, 



P{|G,,(z)-EGii {z)\ >5} <exp 



\B 



N\ 



(18) 



Proof: By proceeding in the same fashion as in the proof of LemmalU we compute 

dC,, 



dB 



GixGyi . 



xy 



Since G is symmetric, this matrix has rank 1 and its only eigenvalue equals {Gai)^ 
Therefore, 



dGi 



dB. 



xy 



< 



IGeJp < IIGII^ < ^ 



,2 ■ 



Hence, we can estimate 



WdxG. 



< 



2\\B 



N\ 



It 2 — 2 



< 



The rest of the proof is similar to the proof of Lemma IH □ 

Lemma 18. Let z = E+ir], 1 » S> N^^/^ , and assume that(i) max {dL{fiAN^ l^a),dL{fj,BM^ l^l^)} 
N^'^^^, (ii) the pair of probability measures (^ail^p) is well-behaved at E, and (Hi) max {\\ A j<4\\ , ||i?7v||} < 
K for all N . Then there exists C > such that for all 6 > 0: 

P{\Gii{z)\>G + 5}< exp (-cJ^iV . 
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Proof: This is a direct consequence of Lemmas [T6l and [TTl (In Lemma [TT] we use rj = 

5. Concluding Remarks 

We showed that there exists a function Sb,n (z) that have the following properties: 

(i) Sb,n (z) is analytical in C+; 

(ii) If Imz > cN~^/^ , and N is sufficiently large then 

Moreover, if jiA^ = /"Biv,! "^^en IvclSb^n {z) > for all z with Imz > 0, and therefore 
5'_B,ivmaps to C+; 

(iii) if Ra,n (z) is defined as 

Ra,n (z) = Ghn (z) - Gan (z + Sb,n (z)) , 

then \\Ra,n iz)\\ ^ provided that Imz > cN''^/^. 

One remaining question is whether it is true that Sb,n (z) maps C+ to C+ without the 
additional assumption that /^Ajv = /^Sjv,- 

By using this result, we showed that a local limit law holds for the empirical eigenvalue 
measure with the window length N~^^^. This result improves over the result in Q. 
However, it is still far form the expected optimal result with the window length A^^^+^. 

Finally, we used the subordination property to show that the localization length of eigen- 
vectors is greater than N^/^. Again, this result is not optimal and the expected localization 
length is O {N) (complete delocalization). 
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